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Capítulo 1 


Empezando con Fortran 

1.1. Observaciones 


Fortran es un lenguaje utilizado ampliamente en la comunidad científica debido a su idoneidad 
para el cálculo numérico. Particularmente atractiva es su intuitiva notación matricial, que facilita 
la escritura de cómputos vectorizados rápidos. 

A pesar de su edad, Fortran todavía está activamente desarrollado, con numerosas implementa- 
ciones, incluyendo GNU, Intel, PGI y Gray. 


1.2. Versiones 


Versión 
FORTRAN66 
FORTRAN?? 
Fortran 90 
Fortran 95 
Fortran 2003 
Fortran 2008 


Nota 

Primera estandarización por ASA (ahora ansí) 
Forma fija, histórica 

Forma libre, norma ISO, operaciones matriciales. 
Procedimientos puros y elementales 
Programación orientada a objetos 
Go-matrices 


Lanzamiento 

?-03-1996 

15-04-19?8 

15-04-1991 

15-6-199? 

4-04-2004 

10-09-2010 


1.3. Instalación o configuración 


Fortran es un lenguaje que puede compilarse utilizando compiladores suministrados por muchos 
proveedores. Diferentes compiladores están disponibles para diferentes plataformas de hardware y 
sistemas operativos. Algunos compiladores son software libre, algunos se pueden utilizar de forma 
gratuita y otros requieren la compra de una licencia. 

El compilador gratuito más común de Fortran es GNU Fortran o gfortran. El código fuente 
está disponible en GNU como parte de GGG, la colección del compilador GNU. Los binarios para 
muchos sistemas operativos están disponibles en https://gcc.gnu.org/wiki/GEortranBinaries . Las 
distribuciones de Linux a menudo contienen gfortran en su gestor de paquetes. 

Otros compiladores están disponibles por ejemplo: 
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EKOPath por PathScale 
LLVM (backend vía DragonEgg) 

Oracle Developer Studio 

Absoft Eortran Compiler 
Compilador Intel Eortran 
Compilador NAG Eortran 
Compiladores de IGP 

En los sistemas HPC, a menudo hay compiladores especializados disponibles por el proveedor 
del sistema como, por ejemplo, los compiladores IBM o Cray . 

Todos estos compiladores soportan el estándar Eortran 95. ACM Eortran Eorum ofrece una des¬ 
cripción general sobre el estado de Eortran 2003 y el estado de Eortran 2008 de varios compiladores 
y está disponible en la Wiki de Eortran. 

1.4. Hola Mundo 

Cualquier programa de Eortran tiene que incluir el end como última declaración. Por lo tanto, 
el programa Eortran más simple se ve así: 

end 

Aquí hay algunos ejemplos de programas ”hola mundo”: 

print *, "Helio, world" 
end 

Con declaración de write 

writeC*,*) "Helio, word" 
end 


Para mayor claridad, ahora es común usar la declaración de program para iniciar un programa 
y darle un nombre. La declaración end puede referirse a este nombre para que sea obvio a que 
se refiere, y permitir que el compilador verifique que el código sea correcto. Además, todos los 
programas de Eortran deben incluir una declaración implicit none. Por lo tanto, un programa 
mínimo de Eortran debería tener el siguiente aspecto: 

program helio 
implicit none 
writeC*,*) ’Hello world!’ 
end program helio 

El siguiente paso lógico a partir de este punto es cómo ver el resultado del programa helio 
world. Esta sección muestra cómo lograr eso en un entorno similar a Linux. Suponemos que tiene 
algunas nociones básicas de los comandos de shell, principalmente sabe cómo llegar al terminal de 
shell. También asumimos que ya ha configurado su entorno f ortran. Utilizando su editor de texto 
preferido (notepad, notepad-|--|-, vi, vim, emacs, gedit, kate, etc.), guarde el programa de saludo de 
arriba (copie y pegue) en un archivo llamado helio.f90 en su directorio de inicio, helio.f90 es 
su archivo fuente. Luego vaya a la línea de comandos y navegue hasta el directorio (¿directorio de 
inicio?) Donde guardó su archivo fuente, luego escriba el siguiente comando: 

>gfortran -o helio helio.f90 
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Acabas de crear tu programa ejecutable helio world. En términos técnicos, acaba de compilar su 
programa. Para ejecutarlo, escriba el siguiente comando: 

>./helio 

Debería ver la siguiente línea impresa en su terminal shell. 

> Helio world! 

Enhorabuena, acaba de escribir, compilar y ejecutar el programa ”Hello World”. 

1.5. Ecuación cuadrática 

Hoy en día Eortran se utiliza principalmente para el cálculo numérico. Este ejemplo muy simple 
ilustra la estructura básica del programa para resolver ecuaciones cuadráticas: 

program quadratic 
! a coimnent 

Ishould be present in every sepárate program unit 
implicit none 

real :: a, b, c 
real :: discriminant 
real :: xl, x2 

print *, "Enter the quadratic equation coefficients a, b and c:" 
read *, a, b, c 

discriminant = b**2 - 4*a*c 

if ( discriminant>0 ) then 

xl = ( -b + sqrt(discriminant)) / (2 * a) 
x2 = ( -b - sqrt(discriminant)) / (2 * a) 
print "Real roots:" 
print *, xl, x2 

I Comparison of floating point numbers for equality is often not recommended. 

I Here, it serves the purpose of illustrating the "else if" construct. 
else if ( discriminant==0 ) then 

xl = - b / (2 * a) 
print *, "Real root:" 
print *, xl 
else 

print *, "No real roots." 
end if 

end program quadratic 
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1.6. Insensibilidad al caso 

Las letras mayúsculas y minúsculas del alfabeto son equivalentes en el conjunto de caracteres 
de Fortran. En otras palabras, Fortran no distingue entre mayúsculas y minúsculas . Este compor¬ 
tamiento está en contraste con los lenguajes que distinguen entre mayúsculas y minúsculas, como 
C y muchos otros. 

Como consecuencia, las variables a y A son la misma variable. En principio se podría escribir un 
programa de la siguiente manera 

pROgrAm MYproGRaM 

enD mYPrOgrAM 
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Capítulo 2 


Alternativas modernas a las 
características históricas. 

2.1. Tipos de variables implícitas 

Cuando Fortran se desarrolló originalmente, la memoria era muy importante. Las variables y los 
nombres de los procedimientos podrían tener un máximo de 6 caracteres, y las variables a menudo se 
escribían de forma implícita. Esto significa que la primera letra del nombre de la variable determina 
su tipo. 

variables que comienzan con i, j , . . . ,n son integer 
todo lo demás (a,b, . . . ,h,yo,p,. . . ,z) son real 

Programas como los siguientes son aceptables Fortran: 

program badbadnotgood 

j = 4 

key = 5 ! only the first letter determines the type 
X = 3.142 

print*, "j = ", j, "key = ", key, "x = ", x 
end program badbadnotgood 

Incluso puede definir sus propias reglas implícitas con la declaración implicit: 

! all variables are real by default 
implicit real (a-z) 

o 


! variables starting with x, y, z are complex 

! variables starting with c, s are character with length of 4 bytes 
! and all other letters have their default implicit type 
implicit complex (x,y,z), character*4 (c,s) 

La tipificación implícita ya no se considera la mejor práctica. Es muy fácil cometer un 
error al utilizar la escritura implícita, ya que los errores tipográficos pueden pasar desapercibidos, 
por ejemplo. 




program oops 

real :: somelongandcomplicatedname 


cali expensive_subroutine(somelongandcomplEcatedname) 
end program oops 

Este programa estará felizmente ejecutado y hará lo incorrecto. 

Para desactivar la escritura implícita, se puede usar la instrucción implicit none. 

program much_better 

implicit none 

integer :: j = 4 

real :: x = 3.142 

print*, "j = ", j, "x = ", X 

end program much_better 

Si hubiéramos usado implicit none en el programa oops arriba, el compilador se habría dado 
cuenta de inmediato y hubiese producido un error. 


2.2. Declaración IF aritmética 

La instrucción aritmética if permite usar tres ramas dependiendo del resultado de una expresión 
aritmética 

if (arith_expr) labell, label2, labelS 

Esta instrucción if transfiere el control de flujo a una de las etiquetas de un código. Si el resultado 
de arith_expr es negativo, labell está involucrado, si el resultado es cero, se usa label2 y si el 
resultado es positivo, se aplica la última labelS. La IF Aritmética requiere las tres etiquetas, pero 
permite la reutilización de las etiquetas, por lo tanto, esta declaración se puede simplificar a una 
rama dos if. 

Ejemplos: 

if (N * N - N / 2) 130, 140, 130 


if (X) 100, lio, 120 


Ahora esta función es obsoleta con la misma funcionabilidad ofrecida por la declaración if y la 
construcción if-else. Por ejemplo, el fragmento: 


if (X) 100, 
100 print*, 
goto 200 
lio print*, 
goto 200 
120 print*, 
200 continué 


lio, 120 

"Negative" 

"Zero" 

"Positive" 


puede ser escrito como la construcción if-else 
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if (X<0) then 
print*, "Negative" 
else if (X==0) then 
print*, "Zero" 
else 

print*, "Positive" 
end if 

Un sustituto if de 

if (X) 100, 100, 200 

100 print *, "Negative or zero" 

200 continué 

tal vez 

if (X<=0) print*, "Negative or zero" 

2.3. Construcción de DO no en bloque 

La construcción do de no bloque parece como 

integer i 
do 100, i=l, 5 
100 print *, i 

Es decir, donde la instrucción de terminación etiquetada no es una instrucción de continué. 
Hay varias restricciones en la declaración que pueden usarse como la declaración de terminación y 
todo es generalmente muy confuso. 

Dicha construcción no de bloque se puede reescribir en forma de bloque como 

integer i 
do 100 i=l,5 
print *, i 
100 continué 

o mejor, usando una sentencia de end do en la terminación, 

integer i 
do i=l,5 
print *, i 
end do 

2.4. Retorno alternativo 

El retorno alternativo es una facilidad para controlar el flujo de ejecución en el retorno desde 
una subrutina. A menudo se utiliza como una forma de manejo de errores: 

real x 

cali sub(x, 1, *100, *200) 
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print*, "Success:", x 
stop 

100 print*, "Negativo input valué" 
stop 

200 print*, "Input valué too large" 
stop 

end 

subroutine sub(x, i, *, *) 
real, intent(out) :: x 
integer, intent(in) :: i 
if (i<0) return 1 
if (i>10) return 2 

X = i 

end subroutine 

El retorno alternativo está marcado por los argumentos * en la lista de argumentos ficticios de 
subrutina. 

En la declaración de cali anterior *100 y *200 refieren a las declaraciones etiquetadas 100 y 
200 respectivamente. 

En la subrutina, las declaraciones de return correspondientes a la devolución alternativa tienen 
un número. Este número no es un valor de retorno, pero denota la etiqueta proporcionada a la 
cual se pasa la ejecución en el retorno. En este caso, el returnul pasa la ejecución a la declaración 
etiquetada 100 y el returnu2 ejecución pasa a la declaración etiquetada 200. Una declaración de 
return sin adornos, o la finalización de la ejecución de subrutina sin una declaración de return, la 
ejecución de aprovado inmediatamente después de la instrucción de llamada. 

La sintaxis alternativa de retorno es muy diferente de otras formas de asociación de argumentos 
y la instalación introduce un control de flujo contrario a los gustos modernos. El control de flujo 
más agradable se puede administrar con la devolución de un código de "estado” entero. 

real x 

integer status 

cali sub(x, 1, status) 
select case (status) 
case (0) 

print*, "Success:", x 
case (1) 

print*, "Negative input valué" 
case (2) 

print*, "Input valué too large" 
end select 

end 

subroutine sub(x, i, status) 
real, intent(out) :: x 
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integer, intent(in) :: i 
integer, intent(out) :: status 

status = 0 

if (i<0) then 
status = 1 
else if (i>10) 
status = 2 
else 
X = i 
end if 

end subroutine 

2.5. Forma de fuente fija 

Fortran originalmente fue diseñado para un formato de formato fijo basado en una tarjeta 
perforada de 80 columnas: Estos fueron creados en una máquina perforadora de tarjetas, muy 



Figura 2.1: Tarjeta perforada 


similar a esto: 



Figura 2.2: 
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El formato, como se muestra en la tarjeta de muestra ilustrada, tenía las primeras cinco columnas 
reservadas para las etiquetas de estados de cuenta. La primera columna se usó para denotar los 
comentarios mediante una letra C. La sexta columna se usó para denotar una continuación de la 
declaración (insertando cualquier carácter que no sea un cero ’O’). Las últimas 8 columnas se usaron 
para la identificación y secuenciación de las tarjetas, ¡lo cual fue muy valioso si dejaste caer el mazo 
de cartas en el suelo! La codificación de caracteres para tarjetas perforadas tenía solo un conjunto 
limitado de caracteres y solo estaba en mayúsculas. Como resultado, los programas de Fortran se 
parecían a esto: 


DIMENSION A(10) 00000001 

C TRIS IS A COMMENT STATEMENT TO EXPLAIN TRIS EXAMPLE PROGRAM 00000002 

WRITE (6,100) 00000003 

100 FORMA!(169RTRIS IS A RATRER LONG STRING BEING OUTPUT WRICR GOES 0VE00000004 
IR MORE IRAN ONE LINE, AND USES IRE STATEMENT CONTINUATION MARKER IN00000005 
2C0LUMN 6, AND ALSO USES ROLLERITR STRING FORMA!) 00000006 

STOP 00000007 

END 00000008 

El carácter de espacio también se ignoró en todas partes, excepto dentro de una constante de carácter 


Hollerith (como se muestra arriba). Esto significaba que los espacios podían aparecer dentro de 
palabras y constantes reservadas, o se podían perder por completo. Esto tuvo el efecto secundario 
de algunas afirmaciones bastante engañosas como: 

DO 1 I = 1.0 

es una asignación a la variable DO II mientras que: 


DOll = 1,0 

es en realidad un bucle de DO en la variable I 

Fortran moderno ahora no requiere esta forma fija de entrada y permite la forma libre usando 
cualquier columna. Los comentarios ahora están indicados por un ! que también se puede añadir a 
una línea de declaración. Los espacios ahora no están permitidos en ningún lugar y deben usarse 
como separadores, como en la mayoría de los otros idiomas. El programa anterior se podría escribir 
en Fortran moderno como: 


! This is a coimnent statement to explain this example program 
Print *,"TRIS IS A RATRER LONG STRING BEING OUTPUT WRICR no longer 
GOES OVER MORE IRAN ONE LINE, AND does not need to USE IRE STATEMENT 
CONTINUATION MARKER IN COLUMN 6, or the ROLLERITR STRING FORMA!" 

Aunque la continuación de estilo antiguo ya no se usa, el ejemplo anterior ilustra que todavía se 
producirán declaraciones muy largas. Modern Fortran usa un símbolo & al final y al comienzo de la 
continuación. Por ejemplo, podríamos escribir lo anterior en una forma más legible: 

! This is a comment statement to explain this example program 
Print *,"TRIS IS A RATRER LONG STRING BEING OUTPUT WRICR still & 

&G0ES OVER MORE IRAN ONE LINE, AND does need to USE IRE STATEMENT & 

&C0NTINUATI0N notation" 
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2.6. Bloques comunes 

En las primeras formas de Fortran, el único mecanismo para crear un almacén de variables global 
visible desde subrutinas y funciones es usar el mecanismo de bloque COMMON . Esto permitió que las 
secuencias de variables fueran nombres y se compartieran en común. 

Además de los bloques comunes con nombre, también puede haber un bloque común en blanco 
(sin nombre). 

Un bloque común en blanco podría ser declarado como 
common i, j 

mientras que las variables bloque nombradas podrían ser declaradas como 
common /variables/ i, j 

Como ejemplo completo, podríamos imaginar una tienda de almacenamiento dinámico utilizada por 
rutinas que pueden agregar y eliminar valores: 

PROGRÁM STACKING 

COMMON /HEAP/ ICOUNT, ISTACK(1023) 

ICOUNT = 0 
READ IVAL 
CALL PUSH(IVAL) 

CALL POP(IVAL) 

END 

SUBROUTINE PUSH(IVAL) 

COMMON /HEAP/ ICOUNT, ISTACK(1023) 

ICOUNT = ICOUNT + 1 
ISTACKCICOUNT) = IVAL 
RETURN 
END 

SUBROUTINE POP(IVAL) 

COMMON /HEAP/ ICOUNT, ISTACK(1023) 

IVAL = ISTACK(ICOUNT) 

ICOUNT = ICOUNT - 1 

RETURN 

END 

Las declaraciones comunes se pueden usar para declarar implícitamente el tipo de una variable y para 
especificar el atributo de dimensión. Este comportamiento solo es a menudo una fuente suficiente 
de confusión. Además, la asociación de almacenamiento implícita y los requisitos para definiciones 
repetidas en las unidades del programa hacen que el uso de bloques comunes sea propenso a errores. 

Finalmente, los bloques comunes están muy restringidos en los objetos que contienen. Por ejem¬ 
plo, una matriz en un bloque común debe ser de tamaño explícito; objetos asignables no pueden 
ocurrir; Los tipos derivados no deben tener inicialización por defecto. 

En Fortran moderno, este uso compartido de variables se puede manejar mediante el uso de 
módulos. El ejemplo anterior se puede escribir como: 

module heap 
implicit none 
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! In Fortran 2008 all module variables are implicitly saved 
integer, save :: count = 0 
integer, save :: stack(1023) 
end module heap 

program stacking 
implicit none 
integer val 
read *, val 
cali push(val) 
cali pop(val) 

contains 

subroutine push(val) 

use heap, only : count, stack 

integer val 

count = count + 1 

stack(count) = val 

end subroutine push 

subroutine pop(val) 

use heap, only : count, stack 

integer val 

val = stack(count) 

count = count - 1 

end subroutine pop 

end program stacking 

Los bloques comunes nombrados y en blanco tienen comportamientos ligeramente diferentes. 

Nota: 

■ Los objetos en bloques comunes nombrados pueden definirse inicialmente; 

■ Los objetos en blanco no serán comunes. 

■ Los objetos en bloques comunes en blanco se comportan como si el bloque común tuviera el 
atributo de save; 

■ los objetos en bloques comunes nombrados sin el atributo de save pueden volverse indefinidos 
cuando el bloque no está en el alcance de una unidad de programa activa 

Este último punto puede contrastarse con el comportamiento de las variables del módulo en el 
código moderno. Todas las variables de módulo en Fortran 2008 se guardan implícitamente y no se 
vuelven indefinidas cuando el módulo queda fuera del alcance. Antes de que las variables del módulo 
Fortran 2008, como las variables en bloques comunes con nombre, también se volvian indefinidas 
cuando el módulo quedaba fuera del alcance. 

2.7. GOTO asignado 

El GOTO asignado utiliza una variable entera a la que se asigna una etiqueta de declaración 
utilizando la instrucción ASSIGN. 
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100 CONTINUE 


ASSIGN 100 TO ILABEL 


GOTO ILABEL 

El GOTO asignado es obsolescente en Fortran 90 y se elimina en Fortran 95 y posteriores. Se 
puede evitar en el código moderno utilizando procedimientos, procedimientos internos, indicadores 
de procedimientos y otras características. 

2.8. GOTO computado 

El GOTO computado permite la bifurcación del programa según el valor de una expresión entera. 
GOTO (label_l, label_2,... label_n) scalar-integer-expression 

Si scalar-integer-expression es igual a 1, el programa continúa en la etiqueta label_l , si es 
igual a 2, va a label_2 y así sucesivamente. Si es menos de 1 o mayor que n programa continúa en 
la siguiente línea. 

Ejemplo: 

ivar = 2 


GOTO (10, 20, 30, 40) ivar 
saltará a la etiqueta de declaración 20. 

Esta forma de goto es obsolescente en Fortran 95 y versiones posteriores, siendo reemplazada 
por la construcción selectucase 

2.9. Especificaciones de formato asignado 

Antes de Fortran 95 era posible usar formatos asignados para entrada o salida. Gonsiderar 

integer i, fmt 
read *, i 

assign 100 to fmt 
if (i<100000) assign 200 to fmt 

print fmt, i 

100 format (’’This is a big number’’, 110) 

200 format (’’This is a small number’’, 16) 

end 
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La declaración de assign asigna una etiqueta de declaración a una variable entera. Esta variable 
entera se usa más tarde como el especificador de formato en la declaración de print. 

Dicha asignación de especificador de formato se eliminó en Fortran 95. En su lugar, un código 
más moderno puede usar alguna otra forma de control de flujo de ejecución 

integer i 
read i 

if (i<100000) then 
print 100, i 
else 

print 200, i 
end if 

100 format ("This is a big number", 110) 

200 format ("This is a small number", 16) 

end 

o se puede usar una variable de carácter como el especificador de formato 

character(29), target :: big_fmt=’("This is a big number", 110)’ 
characterOO), target :: small_fmt=’("This is a small number", 16)’ 
character(:), pointer :: fmt 

integer i 
read *, i 

fmt=>big_fmt 

if (i<100000) fmt=>small_fmt 

print fmt, i 

end 

2.10. Sentencias de Funciones 

Considerar el programa 

implicit none 
integer f, i 
f(i)=i 


print *, f(l) 
end 

Aquí f es una función de declaración. Tiene un tipo de resultado de entero, tomando un argumento 
ficticio de entero. Q 

’Eii los ejemplos de código antiguo anterior, no sería extraño ver que los argumentos ficticios de una función de 
declaración se escriben de forma implícita, incluso si el resultado tiene un tipo explícito. 
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Dicha función de declaración existe dentro del alcance en el que se define. En particular, tiene 
acceso a variables y constantes con nombre accesibles en ese ámbito. 

Sin embargo, las funciones de declaración están sujetas a muchas restricciones y son potencial¬ 
mente confusas (mirando a simple vista como una instrucción de asignación de elementos de matriz). 
Las restricciones importantes son: 

■ El resultado de la función y los argumentos ficticios deben ser escalares. 

■ Los argumentos ficticios están en el mismo ámbito que la función. 

■ Eunciones de declaración no tienen variables locales. 

■ Las funciones de sentencias no se pueden pasar como argumentos reales. 

Los beneficios principales de las funciones de sentencias son repetidos por las funciones internas, 
implicit none 

print *, f(l) 

contains 

integer function f(i) 
integer i 
f = i 

end function 
end 

Las funciones internas no están sujetas a las restricciones mencionadas anteriormente, aunque qui¬ 
zás vale la pena señalar que un subprograma interno puede no contener un subprograma interno 
adicional (pero puede contener una función de declaración). 

Las funciones internas tienen su propio alcance, pero también tienen una asociación de host 
disponible. 
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Capítulo 3 

Matrices 


3.1. Notación básica 

Cualquier tipo se puede declarar como una matriz utilizando el atributo de dimensión o simple¬ 
mente indicando directamente la dimensión o dimensión de la matriz: 

! One dimensional array with 4 elementa 
integer, dimension(4) :: foo 

! Two dimensional array with 4 rows and 2 columna 
real, dimension(4, 2) :: bar 

! Three dimensional array 

type(mytype), dimension(6, 7, 8) :: myarray 

! Same as above without using the dimensión keyword 

integer :: foo2(4) 

real :: bar2(4, 2) 

type(mytype) :: myarray2(6, 7, 8) 

La última forma de declarar una matriz multidimensional, permite la declaración de matrices de 
diferente rango / dimensiones del mismo tipo en una línea, como sigue 

real :: pencil(5), píate(3,-2:4), cuboid(0:3,-10:5,6) 

El rango máximo permitido (número de dimensiones) es 15 en el estándar Fortran 2008 y fue 7 
antes. 

Fortran almacena matrices en orden de columnas mayores. Es decir, los elementos de bar se 
almacenan en la memoria de la siguiente manera: 

bar(1,1),bar(2,1),bar(3,1),bar(4,1),bar(1,2),bar(2,2),... 

En Fortran, la numeración de matrices comienza en 1 por defecto, a diferencia de C que comienza 
en 0. De hecho, en Fortran, puede especificar explícitamente los límites superior e inferior de cada 
dimensión: 

integer, dimensión(7:12, -3:-l) :: geese 
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Esto declara una matriz de forma (6,3), cuyo primer elemento son los geese(7,-3). 

Los límites inferior y superior a lo largo de las 2 (o más) dimensiones pueden accederse mediante 
las funciones intrínsecas ubound y Ibound. De hecho, IboundCgeese ,2) devolvería -3, mientras que 
uboundCgeese, 1) devolvería 12. 

Se puede acceder al tamaño de una matriz mediante el size función intrínseca. Por ejemplo, 
size(geese,udiinu=ul) devuelve el tamaño de la primera dimensión que es 6. 

3.2. Arreglos asignables 

Las matrices pueden tener el atributo asignable: 

! One dimensional allocatable array 
integer, dimension(:), allocatable :: foo 
! Two dimensional allocatable array 
real, dimension(:,:), allocatable :: bar 

Esto declara la variable pero no le asigna ningún espacio. 

! We can specify the bounds as usual 
allocate(foo(3:5)) 

! It is an error to allocate an array twice 
! so check it has not been allocated first 
if (.not. allocated(foo)) then 
allocate(bar(10, 2)) 
end if 

Una vez que una variable ya no es necesaria, se puede desasignar: 
deallocate(foo) 

Si por alguna razón falla una instrucción de allocate, el programa se detendrá. Esto se puede 
evitar si el estado se comprueba mediante la palabra clave stat: 

real, dimension(:), allocatable :: geese 
integer :: status 

allocate(geese(17), stat=status) 
if (stat /= 0) then 

print*, "Something went wrong trying to allocate ’geese’" 
stop 1 
end if 

El deallocate declaración tiene stat de palabras clave también: 
deallocate (geese, stat=status) 

status es una variable entera cuyo valor es 0 si la asignación o desasignación fue exitosa. 
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3.3. Constructores de matrices 


Se puede crear un valor de matriz de rango 1 utilizando un constructor de matriz, con la sintaxis 

(/ ... /) 

[ ... ] 

La forma [. . .] se introdujo en Fortran 2003 y generalmente es considerado como más fáciles de 
interpretar, sobre todo en expresiones complejas. Este formulario se utiliza exclusivamente en este 
ejemplo. Los valores que aparecen en un constructor de matriz pueden ser valores escalares, valores 
de matriz o bucles implícitos. 

Los parámetros de tipo y tipo de la matriz construida coinciden con los de los valores en el 
constructor de matriz 


[1, 2, 3] 
[1., 2., 3.] 
["A", "B"] 


A rank-1 length-3 array of default integer type 

A rank-1 length-3 array of default real type 

A rank-1 length-2 array of default character type 


integer, parameter :: A = [2, 4] 

[1, A, 3] ! A rank-1 length-4 array of default integer type, with A’s elemente 


integer i 

[1, (i, i=2, 5), 6] ! A rank-1 length-6 array of default integer type with an implied-do 


En los formularios anteriores, todos los valores dados deben ser del mismo tipo y tipo de parámetro. 
Los tipos de mezcla, o parámetros de tipo, no están permitidos. Los siguientes ejemplos no son 

válidos. 

[1,2.] ! INVALID: Mixing integer and default real 

[leO, 2d0] ! INVALID: Mixing default real and double precisión 

[1., 2._dp] ! INVALID: Allowed only if kind ‘dp‘ corresponda to default real 

["Helio", "Frederick"] ! INVALID: Different length parameters 

Para construir una matriz utilizando diferentes tipos, se dará una especificación de tipo para la 
matriz 


[integer :: 1, 2., 3d0] ! A default integer array 

[real(dp) :: 1, 2, 3._sp] ! A real(dp) array 

[character(len=9)::"Helio", "Frederick"] ! A length-2 array of length-9 characters 

Esta última forma para matrices de caracteres es especialmente conveniente para evitar el espacio 
de relleno, como la alternativa 

["Helio ", "Frederick"] ! A length-2 array of length-9 characters 

El tamaño de una matriz denominada constante puede estar implícito en el constructor de la matriz 
que se utiliza para establecer su valor 

integer, parameter :: ids(*) = [1, 2, 3, 4] 

y para tipos parametrizados por longitud, el parámetro de longitud puede ser asumido 
character(len=*), parameter :: namesC*) = [character(3) :: "Me", "You", "Her"] 

La especificación de tipo también se requiere en la construcción de matrices de longitud cero. Desde 


21 






[ ] ! Not a valid array constructor 

los parámetros de tipo y tipo no se pueden determinar a partir del conjunto de valores no existentes. 
Para crear una matriz de enteros por defecto de longitud cero: 

[integer :: ] 

Los constructores de matrices construyen solo matrices de rango 1. A veces, como al establecer el 
valor de una constante nombrada, también se requieren arreglos de rango más alto en una expresión. 
Las matrices de rango más alto se pueden tomar del resultado de la reshape con una matriz de rango 
1 construida 

integer, parameter :: multi_rank_ids(2,2) = RESHAPE([1,2,3,4], shape=[2,2]) 

En un constructor de matriz, los valores de la matriz en orden de elementos con cualquier matriz en 
la lista de valores es como si los elementos individuales se dieran en orden de elementos de matriz. 
Así, el ejemplo anterior 

integer, parameter :: A = [2, 4] 

[1, A, 3] ! A rank-1 length-4 array of default integer type, with A’s elemente 

es equivalente a 

[1, 2, 4, 3] ! With the array written out in array element order 

En general, los valores en el constructor pueden ser expresiones arbitrarias, incluidos los constructo¬ 
res de matrices anidadas. Para que un constructor de matriz de este tipo cumpla ciertas condiciones, 
como ser una expresión constante o de especificación, las restricciones se aplican a los valores cons¬ 
tituyentes. 

Aunque no es un constructor de matriz, ciertos valores de matriz también pueden crearse con¬ 
venientemente usando la función intrínseca de spread . Por ejemplo 

[(0, i=l,10)] ! An array with 10 default integers each of valué 0 

También es el resultado de la referencia de la función. 

SPREAD(0, 1, 10) 

3.4. Especificación de la naturaleza de la matriz: rango y forma 

El atributo de dimensión en un objeto especifica que ese objeto es una matriz. Hay, en Eortran 
2008, cinco naturalezas de matrices^ 

■ forma explícita 

■ forma asumida 

■ tamaño asumido 

■ forma diferida 

■ forma implícita 

^Una Especificación Técnica que extiende Fortran 2008 agrega una sexta naturaleza de matriz: rango asumido. 
Esto no está cubierto aquí. 
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Toma las tres matrices de rango 1 Se pueden escribir de manera equivalente como 
integer a, b, c 

dimensiones) a ! Explicit shape (default lower bound 1), extent 5 

dimensionC:) b ! Assumed or deferred shape 

dimensionC*) c ! Assumed size or implied shape array 

Con estos se puede ver que se requiere un contexto adicional para determinar completamente la 
naturaleza de una matriz. 

3.5. Forma explícita 

Una matriz de forma explícita es siempre la forma de su declaración. A menos que la matriz 
se declare como local a un subprograma o construcción de block, los límites que definen la forma 
deben ser expresiones constantes. En otros casos, una matriz de forma explícita puede ser un objeto 
automático, utilizando extensiones que pueden variar en cada invocación de un subprograma o 
block. 

subroutine sub(n) 
integer, intent(in) :: n 

integer a(5) ! A local explicit shape array with constant bound 

integer b(n) ! A local explicit shape array, automatic object 

end subroutine 
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3.6. Forma asumida 


Una matriz de forma asumida es un argumento ficticio sin el atributo allocatable o pointer. 
Dicha matriz toma su forma del argumento real con el que está asociada. 

integer a(5), b(10) 

cali sub(a) ! In this cali the duiiuny argument is like x(5) 

cali sub(b) ! In this cali the dummy argument is like x(10) 

contains 

subroutine sub(x) 

integer x(:) ! Assumed shape dummy argument 

end subroutine sub 

end 

Cuando un argumento ficticio ha tomado forma, el alcance que hace referencia al procedimiento 
debe tener una interfaz explícita disponible para ese procedimiento. 

3.7. Tamaño asumido 

Una matriz de tamaño asumido es un argumento ficticio que tiene su tamaño asumido a partir 
de su argumento real. 

subroutine sub(x) 

integer x(*) ! Assumed size array 

end subroutine 

Las matrices de tamaños asumidos se comportan de manera muy diferente a las matrices de formas 
asumidas y estas diferencias se documentan en otra parte. 

3.8. Forma diferida 

Una matriz de forma diferida es una matriz que tiene el atributo allocatable o pointer . La 
forma de dicha matriz está determinada por su asignación o asignación de puntero. 

integer, allocatable :: a(:) 
integer, pointer :: b(:) 

3.9. Forma implícita 

Una matriz de formas implícita es una constante con nombre que toma su forma de la expresión 
utilizada para establecer su valor 

integer, parameter :: a(*) = [1,2,3,4] 

Las implicaciones de estas declaraciones matriciales en argumentos ficticios deben documentarse en 
otra parte. 
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1. Una Especificación Técnica que extiende Fortran 2008 agrega una sexta naturaleza de matriz: 
rango asumido. Esto no está cubierto aquí. 

2. Se pueden escribir de manera equivalente como 

integer, dimension(5) :: a 
integer, dimensionC:) :: b 
integer, dimensionC*) :: c 

o 

integer a(5) 
integer b(:) 
integer c(*) 


3.10. Arreglos enteros, elementos de arreglo y secciones de arreglo 

Considere la matriz declarada como 
real x(10) 

Entonces tenemos tres aspectos de interés: 

1. Toda la matriz x; 

2. Elementos de matriz, como x(l) ; 

3. Arreglo de secciones, como x(2:6). 

3.11. Arreglos enteros 

En la mayoría de los casos, toda la matriz x refiere a todos los elementos de la matriz como una so¬ 
la entidad. Puede aparecer en sentencias ejecutables, como printu*5uSUM(x)u,upi'intu*5uSIZE(x) 
o x=l. 

Una matriz completa puede hacer referencia a matrices que no tienen forma explícita (como x 
arriba): 

function f(y) 

real, intent(out) :: y(:) 

real, allocatable :: z(:) 

y = 1. ! Intrinsic assignment for the whole array 

z = [1., 2.,] ! Intrinsic assignment for the whole array, invoking allocation 
end function 

Una matriz de tamaño supuesto también puede aparecer como una matriz completa, pero solo en 
circunstancias limitadas (para ser documentada en otra parte). 
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3.12. Elementos de matriz 


Se hace referencia a un elemento de la matriz para dar índices enteros, uno para cada rango de 
la matriz, que denota la ubicación en la matriz completa: 


real x(5,2) 
x(l,l) = 0.2 
x(2,4) = 0.3 


Un elemento de matriz es un escalar. 


3.13. Secciones de matriz 

Una sección de matriz es una referencia a una serie de elementos (quizás solo uno) de una matriz 
completa, utilizando una sintaxis que involucra dos puntos: 


real x(5,2) 
x(: ,1) = 0. 
x(2,:) = 0. 
x(2:4,l) = 0. 
x( 2:3,1:2) = 0. 
x(l:l,l) = 0. 
x([l,3,5],2) = 0. 


Referring to x(l,l), x(2,l), x(3,l), x(4,l) and x(5,l) 
Referring to x(2,l), x(2,2) 

Referring to x(2,l), x(3,l) and x(4,l) 

Referring to x(2,l), x(3,l), x(2,2) and x(3,2) 
Referring to x(l,l) 

Referring to x(l,2), x(3,2) and x(5,2) 


La forma final anterior utiliza un subíndice vectorial. Esto está sujeto a una serie de restricciones 
más allá de otras secciones de la matriz. 

Cada sección de matriz es en sí misma una matriz, incluso cuando solo se hace referencia a un 
elemento. Eso es x(l:l,l) es una matriz de rango 1 y x(l:l,l:l) es una matriz de rango 2. 

Las secciones de matriz en general no tienen un atributo de toda la matriz. En particular, donde 


real, allocatable :: x(:) 

X = [1,2,3] ! X is allocated as part of the assignment 

X = [1,2,3,4] ! X is dealloacted then allocated to a new shape in the assignment 

la asignación 

x(:) = [1,2,3,4,5] ! This is bad when x isn’t the same shape as the right-hand side 

no está permitido: x(:) , aunque una sección de matriz con todos los elementos de x , no es una 
matriz asignable. 


x(:) = [5,6,7,8] 

está bien cuando x es de la forma del lado derecho. 
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3.14. Componentes de matrices de matrices 


type t 
real y(5) 
end type t 

type(t) x(2) 

También podemos referirnos a matrices completas, elementos de matriz y secciones de matriz en 
configuraciones más complicadas. 

De lo anterior, x es una matriz completa. También tenemos 

x(l)%y ! A whole array 

x(l)%y(l) ! An array element 

x%y(l) ! An array section 

x(l)%y(:) ! An array section 

x([l,2]%y(l) ! An array section 

x(l)%y(l:l) ! An array section 

En tales casos, no se nos permite tener más de una parte de la referencia que consiste en una matriz 
de rango 1. Lo siguiente, por ejemplo, no está permitido 

x%y ! Both the x and y parts are arrays 

x(l: DXy (1:1) ! Recall that each part is still an array section 

3.15. Operaciones de matriz 

Debido a sus objetivos computacionales, las operaciones matemáticas en matrices son sencillas 
en Fortran. 

3.16. Adición y snstracción 

Las operaciones en matrices de la misma forma y tamaño son muy similares al álgebra matricial. 

En lugar de recorrer todos los índices con bucles, se puede escribir suma (y resta): 

real, dimension(2,3) :: A, B, C 
real, diinension(5,6,3) :: D 

A =3. ! Assigning single valué to the whole array 

B =5. ! Equivalent writing for assignment 

C = A + B ! All elementa of C now have valué 8. 

D = A + B ! Compilen will raise an error. The shapes and dimensions are not the same 

Las matrices de corte son también válidas: 
integer :: i, j 

real, dimension(3,2) :: Mat = 0. 

real, dimension(3) :: Ved = 0., Vec2 = 0., Vec3 = 0. 
i = 0 

j = 0 

do i = 1,3 
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do j = 1,2 
Mat(i,j) = i+j 
enddo 
enddo 

Vecl = Mat(:,l) 

Vec2 = Mat(:,2) 

Vec3 = Mat(l:2,l) + Mat(2:3,2) 

3.17. Función 

De la misma manera, la mayoría de las funciones intrínsecas se pueden usar de forma implícita 
asumiendo una operación de componentes (aunque esto no es sistemático): 

real, diinension(2) :: A, B 
A(l) = 6 

A(2) = 44 ! Random valúes 

B = sin(A) ! Identical to B(l) = sin(6), B(2) = sin(44). 

3.18. Multiplicación y división 

Se debe tener cuidado con el producto y la división: las operaciones intrínsecas que usan * y / 
son símbolos en forma elemental: 

real, dimension(2) :: A, B, C 

A(l) = 2 

A(2) = 4 

B(l) = 1 

B(2) = 3 

C = A*B ! Returns C(l) = 2*1 and C(2) = 4*3 

Esto no debe confundirse con operaciones matriciales (ver más abajo). 

3.19. Operaciones matriciales 

Las operaciones matriciales son procedimientos intrínsecos. Por ejemplo, el producto matricial 
de las matrices de la sección anterior se escribe de la siguiente manera: 

real, dimensión(2,1) :: A, B 

real, dimension(l,1) :: C 

A(l) = 2 

A(2) = 4 

B(l) = 1 

B(2) = 3 

C = matmul(transpose(A),B) ! Returns the scalar product of vectors A and B 

Las operaciones complejas permiten encapsular funciones mediante la creación de arreglos tempora¬ 
les. Mientras que algunos compiladores y opciones de compilación lo permiten, esto no se recomienda. 
Por ejemplo, un producto que incluye una transposición matricial puede escribirse: 
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real, dimensión(3,3) :: A, B, C 
A(:) = 4 
B(:) = 5 

C = matmul(transpose(A),matmul(B,matmul(A,transpose(B)))) ! Equivalen! to A~t.B.A.B~T 

3.20. Secciones de matrices avanzadas: tripletes y subíndices vec¬ 
toriales 

Como se mencionó en otro ejemplo , se puede hacer referencia a un subconjunto de los elementos 
de una matriz, llamada sección de matriz. De ese ejemplo podemos tener 

real x(10) 
x(:) = 0. 

x(2:6) = 1. 
x(3:4) = [3., 5.] 

Sin embargo, las secciones de matriz pueden ser más generales que esto. Pueden tomar la forma de 
tripletes de subíndice o subíndices vectoriales. 


3.21. Triplet es de subíndices 

Un subíndice triple toma la forma [boundl]:[bound2][:stride] . Por ejemplo 

real x(10) 
x(l:10) = .. 
x(l:) = ... 
x(:10) = ... 
x(l:6:2) = . 
x(5:l) = ... 
x(5:l:-l) = 
x(::3) = ... 
x(::-3) = . . 


Elemente x(l;, x(2;, ..., x(10; 

The omitted second bound is equivalen! !o the upper, same as above 
The omitted first bound is equivalen! to the lower, same as above 
Elemente x(l), x(3), x(5) 

No elemente: the lower bound is greater than the upper 
Elemente x(5), x(4), x(3), x(2), x(l) 

Elemente x(l), x(4), x(7), x(10), assuming omitted bounds 

No elemente: the bounds are assumed with the first the lower. nesative stride 


Cuando se especifica un paso (que no debe ser cero), la secuencia de elementos comienza con el 
primer límite especificado. Si la zancada es positiva (resp. Negativa), los elementos seleccionados 
siguiendo una secuencia incrementada (resp. Decrementada) por la zancada hasta que el último 
elemento no sea mayor (resp. Más pequeño) que el segundo límite se toma. Si se omite el paso, se 
trata como si fuera uno. 

Si el primer límite es mayor que el segundo límite, y la zancada es positiva, no se especifican 
elementos. Si el primer límite es más pequeño que el segundo límite y la zancada es negativa, no se 
especifican elementos. 

Cabe señalar que x(10:l:-l) no es lo mismo que x(l:10:l) aunque cada elemento de x aparezca 
en ambos casos. 


3.22. Subíndices vectoriales 

Un subíndice vectorial es una matriz de enteros de rango 1. Esto designa una secuencia de 
elementos correspondientes a los valores de la matriz. 
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Elements x(l), x(6), x(4) 
Elements x(2), x(3) and x(4) 
Elements x(2), x(5) and x(2) 


real x(10) 
integer i 
x( [1,6,4]) = ... 
x( [(i,i=2,4)]) = . . . 
print*, x( [2,5,2]) 

Una sección de matriz con un subíndice vectorial está restringida en la forma en que se puede usar: 

puede que no sea un argumento asociado con un argumento ficticio que se define en el proce¬ 
dimiento; puede que no sea el objetivo en una instrucción de asignación de puntero; puede que no 
sea un archivo interno en una declaración de transferencia de datos. Además, tal sección de ma¬ 
triz puede no aparecer en una declaración que implique su definición cuando el mismo elemento se 
selecciona dos veces. Desde arriba: 

print*, x( [2,5,2]) ! Elements x(2), x(5) and x(2) are printed 

x( [2,5,2]) = 1. ! Not permitted: x(2) appears twice in this definition 

3.23. Secciones de matriz de rango más alto 

real x(5,2) 

print*, x(::2,2:1:-1) ! Elements x(l,2), x(3,2), x(5,2), x(l,l), x(3,l), x(5,l) 
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Capítulo 4 

Control de ejecución 


4.1. Construcción IF 

El constructo if (llamado una instrucción IF de bloque en FORTRAN 77) es común en muchos 
lenguajes de programación. Se ejecuta condicionalmente un bloque de código cuando una expresión 
lógica se evalúa como verdadera. 

[ñame:] IF (expr) THEN 
block 

[ELSE IF (expr) THEN [ñame] 
block] 

[ELSE [ñame] 
block] 

END IF [ñame] 
dónde 

nombre - el nombre de la construcción if (opcional) 
expr - una expresión lógica escalar entre paréntesis 
bloque - una secuencia de cero o más declaraciones o construcciones 

Un nombre de construcción al comienzo de una instrucción if then debe tener el mismo valor que 
el nombre de construcción al end if instrucción end if , y debe ser único para la unidad de alcance 
actual. 

En las sentencias if , las ecuaciones (in) y las expresiones lógicas que evalúan una sentencia se 
pueden utilizar con los siguientes operadores: 


LT. 

which is < 

! less than 


LE. 

<= 

! less than or 

equal 

GT. 

> 

! greater than 


GE. 

>= 

! greater than 

or equal 

EQ. 

= 

! equal 


NE. 

/= 

! not equal 


AND. 


! logical and 


OR. 


! logical or 


NOT. 


! negation 
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Ejemplos: 

! simplest form of if construct 
if (a > b) then 
c = b / 2 
end if 

!equivalent example with altérnate syntax 

if(a.gt.b)then 

c=b/2 

endif 

! named if construct 
circle: if (r >= 0) then 
1 = 2 * pi * r 
end if circle 

! complex example with nested if construct 
block: if (a < e) then 
if (abs(c - e) <= d) then 
a = a * c 
else 

a = a * d 
end if 
else 

a = a * e 
end if block 

Un uso histórico de la construcción if encuentra en lo que se denomina una declaración .Aritmética 
if". Sin embargo, dado que esto puede ser reemplazado por construcciones más modernas, no se 
trata aquí. Más detalles se pueden encontrar aquí. 

4.2. Construcción SELECT CASE 

Una construcción de select case condicionalmente ejecuta un bloque de construcciones o decla¬ 
raciones según el valor de una expresión escalar en una declaración de select case . Este constructo 
de control puede considerarse como un reemplazo para el goto computado. 

[ñame:] SELECT CASE (expr) 

[CASE (case-value [, case-value] ...) [ñame] 
block]... 

[CASE DEFAULT [ñame] 
block] 

END SELECT [ñame] 
dónde, 

nombre el nombre de la construcción de selectycase (opcional) 

expr una expresión escalar de tipo entero, lógico o carácter (entre paréntesis) 

bloque una secuencia de cero o más declaraciones o construcciónes 
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Ejemplos: 


! simplest form of select case construct 

select case(i) 

case(:-1) 

s = -1 

case(0) 

s = 0 

case (1:) 

s = 1 

case default 

print "Something strange is happened" 
end select 

En este ejemplo, (:-l) valor de un caso es un rango de valores que coincide con todos los valores 
menores que cero, (0) coincide con los ceros, y (1:) coincide con todos los valores por encima de 
cero, la sección default incluye si otras secciones lo hicieron sin ejecutar. 

4.3. Construcción de bloques DO 

Una construcción do es una construcción de bucle que tiene una serie de iteraciones gobernadas 
por un control de bucle 

integer i 
do i=l, 5 
print *, i 
end do 
print *, i 

En la forma anterior, la variable de bucle i pasa a través del bucle 5 veces, tomando los valores de 
1 a 5 por turno. Una vez que la construcción ha completado, la variable de bucle tiene el valor 6, es 
decir, la variable de bucle se incrementa una vez más después de completar el bucle. 
Más generalmente, la construcción de bucle do se puede entender de la siguiente manera 

integer i, first, last, step 
do i=first, last, step 
end do 

El bucle comienza con i con el valor first, incrementando cada iteración por step hasta que i sea 
mayor que la last (o menor que la last si el tamaño del paso es negativo). 

Es importante tener en cuenta que desde Eortran 95, la variable de bucle y las expresiones de 
control de bucle deben ser enteras. 

Una iteración se puede finalizar de manera prematura con la instrucción de cycle 

do i=l, 5 
if (i==4) cycle 
end do 

y toda la construcción puede cesar la ejecución con la declaración de exit 
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do i=l, 5 
if (i==4) exit 
end do 
print *, i 

Las construcciones do pueden ser nombradas: 

do.name: do i=l, 5 
end do do_name 

Lo cual es particularmente útil cuando hay construcciones anidadas do 

dol: do i=l, 5 
do j=l,6 
if (j==3) cycle 
if (j==4) cycle 
if (i+j==7) cycle dol 
if (i*j==15) exit dol 
end do 
end dol 

Las construcciones do también pueden tener un control de bucle indeterminado, ya sea ”para siem¬ 
pre” o hasta que se cumpla una condición determinada 

integer :: i=0 
do 

i=i+l 

if (i==5) exit 
end do 

o 

integer :: i=0 
do while (i<6) 
i=i+l 
end do 

Esto también permite un bucle do infinito a través de una declaración .true. 

print *,’forever’ 
do while(.true.) 
print =t=,’and ever’ 
end do 

Un constructor do también puede dejar el orden de iteraciones indeterminado. 

do concurrent (i=l:5) 
end do 


This cycles the j construct 
This cycles the j construct 
This cycles the i construct 
This exits the i construct 
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teniendo en cuenta que la forma de control de bucle es la misma que en un control forall . 

Hay varias restricciones en las declaraciones que pueden ejecutarse dentro del rango de una 
construcción do concurrent que están diseñadas para garantizar que no haya dependencias de datos 
entre las iteraciones de la construcción. Esta indicación explícita por parte del programador puede 
permitir una mayor optimización (incluida la paralelización) por parte del compilador, lo que puede 
ser difícil de determinar de otro modo. 

Las variables "privadas"dentro de una interacción se pueden realizar mediante el uso de una 
construcción de block dentro del do concurrent : 

do concurrent (i=l:5, j=2:7) 
block 

real tempval ! This is independent across iterations 
end block 
end do 

Otra forma del bloque do utiliza una instrucción continué etiquetada en lugar de un end do: 

do 100, i=l, 5 
100 continué 

Incluso es posible anidar tales construcciones con una declaración de terminación compartida 

do 100, i=l,5 
do 100, j=l,5 
100 continué 

Ambas formas, y especialmente la segunda (que es obsolescente), generalmente deben evitarse en 
aras de la claridad. 

Finalmente, también es una construcción no-bloque do. Esto también se considera que está 
obsoleto y se describe en otra parte, junto con los métodos para reestructurar como un bloque do. 

4.4. Construcción WHERE 

El where constructo, disponible en Fortran90 representa en adelante un enmascarado do cons¬ 
truir. La declaración de enmascaramiento sigue las mismas reglas de la instrucción if , pero se aplica 
a todos los elementos de la matriz dada. Usando where permite que las operaciones se lleven a 
cabo en una matriz (o múltiples matrices del mismo tamaño), cuyos elementos satisfacen una regla 
determinada. Esto se puede usar para simplificar operaciones simultáneas en varias variables. 
Sintaxis: 

[ñame]: where (mask) 
block 

[elsewhere (mask) 
block] 

[elsewhere 

block] 

end where [ñame] 

Aquí, 

ñame -es el nombre dado al bloque (si es nombrado) 
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mask -es una expresión lógica aplicada a todos los elementos 

block -serie de comando a ejecutar 

Ejemplos: 

! Example variables 
real:: A(5),B(5),C(5) 

A = 0.0 
B = 1.0 

C = [0.0, 4.0, 5.0, 10.0, 0.0] 

! Simple where construct use 
where (C/=0) 

A=B/C 

elsewhere 

A=0.0 

end 

! Named where construct 
Block: where (C/=0) 

A=B/C 

elsewhere 

A=0.0 

end where Block 
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Capítulo 5 


Extensiones de archivo fuente (-f, .f90, 
.f95....) y como se relacionan con el 
compilador. 

5.1. Introducción 

Los archivos de Fortran están bajo una variedad de extensiones y cada uno de ellos tiene un 
significado diferente. Especifican la versión de lanzamiento de Fortran, el estilo de formato de código 
y el uso de directivas de preprocesador similares al lenguaje de programación C. 


5.2. Extensiones y significado 

Las siguientes son algunas de las extensiones comunes que se utilizan en los archivos fuente de 
Fortran y las funcionalidades en las que pueden trabajar. 

F minúscula en la extensión 

Estos archivos no tienen las características de las directivas de preprocesador similares al lenguaje 
de programación C. Se pueden compilar directamente para crear archivos de objetos. Por ejemplo: 
•f, .for, .f95 

F mayúscula en la extensión 

Estos archivos tienen las características de las directivas de preprocesador similares al lenguaje 
de programación C. Los preprocesadores se definen dentro de los archivos o utilizan archivos de 
encabezado como C / C ++ o ambos. Estos archivos deben preprocesarse para obtener los archivos 
de extensión en minúsculas que se pueden usar para compilar. Por ejemplo: .F, .FOR, .F95 

.f, .for, .f77, .ftn 

Se usan para archivos Fortran que usan el formato de estilo Fijo y, por lo tanto, usan la versión 
de lanzamiento de Fortran 77. Como son extensiones en minúsculas, no pueden tener directivas 
de preprocesador. 

.F, .FOR, .F77, .FTN 

Se usan para archivos Fortran que usan el formato de estilo Fijo y, por lo tanto, usan la versión 
de lanzamiento de Fortran 77. Como son extensiones en mayúsculas, pueden tener directivas de 
preprocesador y, por lo tanto, deben procesarse previamente para obtener los archivos de extensión 
en minúsculas. 

.f90, .f95, .f03, .f08 Se usan para archivos Fortran que usan el formato de estilo libre y, por 
lo tanto, usan versiones de Fortran más recientes. Las versiones de lanzamiento están en el nombre. 
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f90 - Fortran 90 
f95 - Fortran 95 
f03 - Fortran 2003 
f08 - Fortran 2008 

Como son extensiones en minúsculas, no pueden tener directivas de preprocesador. 

.F90, .F95, .F03, .F08 Se usan para archivos Fortran que usan el formato de estilo libre 
y, por lo tanto, usan versiones de Fortran más recientes. Las versiones de lanzamiento están en el 
nombre. 

F90 - Fortran 90 
F95 - Fortran 95 
F03 - Fortran 2003 
F08 - Fortran 2008 

Como son extensiones en mayúsculas, tienen directivas de preprocesador y, por lo tanto, deben 
procesarse previamente para obtener los archivos de extensión en minúsculas. 
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Capítulo 6 

i/o 


6.1. Sintaxis 

WRITE(unit num, format num) genera los datos después de los corchetes en una nueva línea. 

READ(unit num, format num) entradas a la variable después de los corchetes 

OPEN(unit num, EILE=file) abre un archivo. (Hay más opciones para abrir archivos, pero no 
son importantes para E / S 

CLOSE(unit num) cierra un archivo. 

6.2. E/S simple 

Como ejemplo de escritura de entrada y salida, tomaremos un valor real y devolveremos el valor 
y su cuadrado hasta que el usuario ingrese un número negativo. 

Como se especifica a continuación, el comando de read toma dos argumentos: el número de 
unidad y el especificador de formato. En el siguiente ejemplo, usamos * para el número de unidad 
(que indica stdin) y * para el formato (que indica el valor predeterminado para reales, en este caso). 
También especificamos el formato para la declaración de print . Alternativamente, se puede usar 
write(*,"The valué....") o simplemente ignorar el formato y tenerlo como 

print *,"The entered valué was ", x," and its square is ",x*x 

lo que probablemente resultará en algunas cadenas y valores espaciados de manera extraña. 

program Simple10 
implicit none 

integer, parameter :: wp = selected_real_kind(15,307) 
real(kind=wp) :: x 

! we’ll loop over until user enters a negativa number 

print ’C'Enter a number >= 0 to see its square. Enter a number < 0 to exit.")’ 
do 

! this reads the input as a double-pricision valué 
read(*,*) x 
if (x < OdO) exit 

! print the entered valué and it’s square 

print ’C'The entered valué was ",fl2.6,", its square is ",f12.6,".’,x,x*x 
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end do 

print ’("Thank you!")’ 
end program SimplelO 

6.3. Leer con un poco de comprobación de errores 

Un ejemplo moderno de Fortran que incluye la verificación de errores y una función para obtener 
un nuevo número de unidad para el archivo. 

module functions 

contains 

function get_new_fileunitO result (f) 
implicit none 

logical :: op 

integer :: f 

f = 1 
do 

inquine(f,opened=op) 
if (op .eqv. .false.) exit 
f = f + 1 
enddo 

end function 
end module 
program file_read 

use functions, only : get_new_fileunit 
implicit none 

integer :: unitno, ierr, readerr 

logical :: exists 

real(kind(0.dO)) :: somevalue 

character(len=128) :: filename 

filename = "somefile.txt" 

inquine(file=trim(filename), exist=exists) 

if (exists) then 

unitno = get_new_fileunitO 

open(unitno, file=trim(filename), action="read", iostat=ierr) 
if (ierr .eq. 0) then 

read(unitno, *, iostat=readerr) somevalue 
if (readerr .eq. 0) then 
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print*, "Valué in file ", trimCfilename), " is ", somevalue 
else 

print*, "Error ", readerr, & 

" attempting to read file ", & 

trim(filename) 

endif 

else 

print*, "Error ", ierr ," attempting to open file ", trim(filename) 

stop 

endif 

else 

print*, "Error -- cannot find file: ", trim(filename) 

stop 

endif 

end program file_read 

6.4. Pasando argumentos de línea de comando 

Donde se admiten los argumentos de la línea de comando, se pueden leer a través del in¬ 
trínseco get_command_argument (introducido en el estándar Fortran 2003). El comando intrínseco 
command_argument_count proporciona una manera de conocer el número de argumentos proporcio¬ 
nados en la línea de comandos. 

Todos los argumentos de la línea de comando se leen como cadenas, por lo que se debe realizar 
una conversión de tipo interna para los datos numéricos. Como ejemplo, este código simple suma 
los dos números proporcionados en la línea de comando: 

PROGRAM cmdlnsum 
IMPLICIT NONE 

CHARACTER(IOO) :: numlchar 
CHARACTER(IOO) :: num2char 


REAE : 

: numl 

REAE : 

: num2 

REAE : 

: numsum 


¡First, make sure the right number of inputs have been provided 
IF(C0MMAND_ARGUMENT_C0UNT().NE.2)TREN 

WRITE(*,*)’ERROR, TWO COMMAND-LINE ARGUMENTE REQUIRED, STOPPING’ 

STOP 

ENDIF 

CALE GET_C0MMAND_ARGUMENT(1,numlchar) Ifirst, read in the two valúes 
CALE GET_C0MMAND_ARGUMENT(2,num2char) 

READ(numlchar,*)numl Ithen, convert them to REAEs 

READ(num2char,*)num2 

numsum=numl+num2 Isum numbers 

WRITE(*,*)numsum Iwrite out valué 
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END PROGRAM 


El argumento del número en get_coinmand_arguinent útilmente entre 0 y el resultado de com- 
mand_argument_count. Si el valor es 0 se proporciona el nombre del comando (si se admite). 

Muchos compiladores también ofrecen intrínsecos no estándar (como getarg ) para acceder a los 
argumentos de la línea de comandos. Como estos no son estándar, se debe consultar la documenta¬ 
ción del compilador correspondiente. 

El uso de get_command_argument puede extenderse más allá del ejemplo anterior con los argu¬ 
mentos de length y status . Por ejemplo, con 

character(5) arg 
integer stat 

cali get_cominand_arguinent (nuinber=l, value=arg, status=stat) 

el valor de stat será -1 si el primer argumento existe y tiene una longitud mayor que 5. Si hay alguna 
otra dificultad para recuperar el argumento, el valor de stat será un número positivo (y arg constará 
completamente de espacios en blanco). De lo contrario su valor será 0 . 

El argumento de length se puede combinar con una variable de carácter de longitud diferida, 
como en el siguiente ejemplo. 

character(:), allocatable :: arg 
integer arglen, stat 

cali get_coimnand_argument(number=l, length=arglen) ! Assume for simplicity success 
allocate (character(arglen) :: arg) 

cali get_cominand_arguinent (number=l, value=arg, status=stat) 
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Capítulo 7 

Interfaces explícitas e implícitas 


7.1. Subprogramas internos/módulos e interfaces explícitas 

Un subprograma (que define un procedimiento ), puede ser una subroutine o una function ; se 
dice que es un subprograma interno si se llama o se invoca desde el mismo program o subprograma 
que lo contains , como sigue 

program my_program 

! declarations 
! executable statements, 

! among which an invocation to 
! internal procedure(s), 
cali my_sub(argl,arg2,...) 
fx = my_fun(xxl,xx2,.. .) 

contains 

subroutine my_sub(al,a2,...) 

! declarations 
! executable statements 
end subroutine my_sub 

function my_fun(xl,x2,...) result(f) 

! declarations 
! executable statements 
end function my_fun 

end program my_program 

En este caso, el compilador sabrá todo sobre cualquier procedimiento interno, ya que trata la unidad 
del programa como un todo. En particular, "verá"la interface del procedimiento, es decir, 

si es una function o subroutine , cuáles son los nombres y propiedades de los argumentos al , a2 , 
xl , x2 , cuáles son las propiedades del resultado f (en el caso de una function ). 

Al conocerse la interfaz, el compilador puede verificar si los argumentos reales ( argl , arg2 , xxl , 
xx2 , fx , ...) pasaron al procedimiento con los argumentos ficticios ( al , a2 , xl , x2 , f , 
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En este caso decimos que la interfaz es explícita . 

Se dice que un subprograma es un subprograma de módulo cuando es invocado por una decla¬ 
ración en el módulo que lo contiene, 

module my_mod 

! declarations 

contains 

subroutine my_mod_sub(bl,b2,...) 

! declarations 
! executable statements 
r = my_mod_fun(bl,b2,...) 
end subroutine my_sub 

function my_mod_fun(yl,y2,...) result(g) 

! declarations 
! executable statements 
end function my_fun 

end module my_mod 

o por una declaración en otra unidad de programa que use s ese módulo, 

program my_prog 

use my_mod 

cali my_mod_sub(...) 

end program my_prog 

Como en la situación anterior, el compilador sabrá todo sobre el subprograma y, por lo tanto, 
decimos que la interfaz es explícita . 


7.2. Subprogramas externos e interfaces implícitas 

Se dice que un subprograma es externo cuando no está contenido en el programa principal, ni 
en un módulo o subprograma adicional. En particular, se puede definir por medio de un lenguaje 
de programación que no sea Eortran. 

Cuando se invoca un subprograma externo, el compilador no puede acceder a su código, por lo 
que toda la información permitida al compilador está implícitamente contenida en la declaración 
del programa que llama y en el tipo una propiedad de los argumentos reales, no los argumentos 
ficticios ( cuya declaración es desconocida para el compilador). En este caso decimos que la interfaz 
es implícita . 

Se puede usar una declaración external para especificar que el nombre de un procedimiento es 
relativo a un procedimiento externo, 

external external_naine_list 
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Pero aun así, la interfaz permanece implícita. 

Se puede utilizar un bloque de interface para especificar la interfaz de un procedimiento 
externo, 

interface 
interface_body 
end interface 

donde interf ace_body es normalmente una copia exacta del encabezado del procedimiento seguido 
de la declaración de todos sus argumentos y, si es una función, del resultado. 

Por ejemplo, para la función WindSpeed 

real function WindSpeedCu, v) 
real, intent(in) :: u, v 
WindSpeed = sqrt(u*u + v*v) 
end function WindSpeed 

Puedes escribir la siguiente interfaz 
interface 

real function WindSpeedCu, v) 
real, intent(in) :: u, v 
end function WindSpeed 
end interface 
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Capítulo 8 

Interoperabilidad C 


8.1. Llamando C desde Fortran 

Fortran 2003 introdujo características de lenguaje que pueden garantizar la interoperabilidad 
entre C y Fortran (y para más idiomas utilizando C como intermediario). A estas funciones se accede 
principalmente a través del módulo intrínseco iso_c_binding: 

use, intrinsic :: iso_c_binding 

La palabra clave intrinsic aquí garantiza que se use el módulo correcto, y no un módulo creado 
por el usuario con el mismo nombre. 

iso_c_binding da acceso a parámetros de tipo de tipo interoperables: 

integer(c_int) :: foo ! equivalent of ’int foo’ in C 

real(c_float) :: bar ! equivalent of ’float bar’ in C 

El uso de parámetros de tipo C garantiza que los datos pueden transferirse entre los programas C 
y Fortran. 

La interoperabilidad de los caracteres C char y Fortran es probablemente un tema por sí mismo 
y, por lo tanto, no se trata aquí. 

Para llamar realmente a una función C desde Fortran, primero se debe declarar la interfaz. Esto 
es esencialmente equivalente al prototipo de la función C, y le permite al compilador conocer el 
número y el tipo de los argumentos, etc. El atributo de bind se usa para decirle al compilador el 
nombre de la función en C, que puede ser diferente al de Eortran. nombre, 
gansos.h 

// Count how many geese are in a given flock 
int howManyGeese(int flock); 

gansos f90 

! Interface to C routine 
interface 

integer(c_int) function how_inany_geese (flock_nuiii) bind(C, ’howManyGeese’) 

! Interface blocks don’t know about their context, 

! so we need to use iso_c_binding to get c_int definition 

use, intrinsic :: iso_c_binding, only : c_int 

integer(c_int) :: flock_num 

end function how_iiiany_geese 

end interface 
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El programa de Fortran debe estar vinculado con la librería C (dependiente del compilador, ¿se 
incluyen aquí?), Que incluye la implementación de howManyGeese() , y luego how_inany_geese() 
puede ser llamado desde Fortran. 


8.2. Estructuras de C en Fortran 

El atributo bind también se puede aplicar a tipos derivados: 
gansos .h 

struct Goose { 
int flock; 
float buoyancy; 

} 

struct Goose goose_c; 
gansos f90 

use, intrinsic :: iso_c_binding, only : c_int, c_float 

type, bind(C) :: goose_t 
integer(c_int) :: flock 
real(c_float) :: buoyancy 
end type goose_t 

type(goose_t) :: goose_f 

Los datos ahora se pueden transferir entre goose_c y goose_f. Las rutinas C que toman argumentos 
de tipo Goose se pueden llamar desde Fortran con type (goose_t). 
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Capítulo 9 

Procedimientos - Funciones y subrutinas. 


9.1. Observaciones 

Las funciones y las subrutinas, junto con los módulos, son las herramientas para dividir un 
programa en unidades. Esto hace que el programa sea más legible y manejable. Cada una de estas 
unidades puede considerarse como parte del código que, idealmente, podría compilarse y probarse 
de forma aislada. Los programas principales pueden llamar (o invocar) a dichos subprogramas 
(funciones o subrutinas) para realizar una tarea. 

Las funciones y las subrutinas son diferentes en el siguiente sentido: 

Las funciones devuelven un solo objeto y, generalmente, no alteran los valores de sus argumentos 
(es decir, ¡actúan como una función matemática); 

Las subroutinas generalmente realizan una tarea más complicada y generalmente alteran sus 
argumentos (si hay alguno presente), así como otras variables (por ejemplo, aquellas declaradas 
en el módulo que contiene la subrutina). 

Las funciones y subrutinas van colectivamente bajo el nombre de procedimientos. (En lo siguiente, 
usaremos el verbo "cali” como sinónimo de "invocar” incluso si técnicamente los procedimientos a 
ser cali son subroutine, mientras que las functions aparecen como parte derecha de la asignación 
o en expresiones). 

9.2. Sintaxis funcional 

Las funciones se pueden escribir utilizando varios tipos de sintaxis. 

function ñame O 

integer ñame 

ñame = 42 

end function 

integer function nameO 

ñame = 42 

end function 

function ñame O result(res) 
integer res 
res = 42 
end function 


48 



Las funciones devuelven valores a través de un resultado de función. A menos que la declaración de 
función tenga una cláusula de result el result la función tendrá el mismo nombre que la función. 
Con result el result la función es el dado por el result. En cada uno de los dos primeros ejemplos 
anteriores, el resultado de la función viene dado por ñame; en el tercero por res. 

El resultado de la función debe definirse durante la ejecución de la función. 

Las funciones permiten utilizar algunos prefijos especiales. 

Eunción pura significa que esta función no tiene ningún efecto secundario: 

puré real function square(x) 
real, intent(in) :: x 
square = x * x 
end function 

La función elemental se define como operador escalar, pero se puede invocar con array como argu¬ 
mento real, en cuyo caso la función se aplicará de forma elemental. A menos que se especifique el 
prefijo impure (introducido en Eortran 2008), una función elemental también es una función pura. 

elemental real function square(x) 
real, intent(in) :: x 
square = x * x 
end function 

9.3. Declaración de devolución 

La declaración de return se puede utilizar para salir de la función y la subrutina. A diferencia 
de muchos otros lenguajes de programación, no se utiliza para establecer el valor de retorno. 


real function f(x) 
real, intent(in) :: x 
integer :: i 

f = X 

do i = 1, 10 

f = sqrt(f) - 1.0 

if (f < 0) then 
f = -1000. 
return 
end if 

end do 

end function 

Esta función realiza un cálculo iterativo. Si el valor de f vuelve negativo, la función devuelve el 
valor -1000. 
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9.4. Procedimientos recursivos 

En Fortran, las funciones y las subrutinas deben declararse explícitamente como recursivas, si 
han de llamarse a sí mismas de nuevo, directa o indirectamente. Por lo tanto, una implementación 
recursiva de la serie Fibonacci podría tener este aspecto: 

recursive function fibonacci(term) result(fibo) 
integer, intent(in) :: term 
integer :: fibo 

if (term <= 1) then 

fibo = 1 

else 

fibo = fibonacci(term-1) + fibonacci(term-2) 
end if 

end function fibonacci 

Otro ejemplo permitido es calcular factorial: 

recursive function factorial(n) result(f) 
integer :: f 

integer, intent(in) :: n 

if(n == 0) then 

f = 1 

else 

f = n * f(n-l) 
end if 

end function factorial 

Para que una función se refiera directamente a sí misma, su definición debe usar el sufijo de result. 
No es posible que una función sea tanto recursive como elemental. 

9.5. La intención de los argumentos ficticios 

El atributo de intent de un argumento ficticio en una subrutina o función declara su uso 
previsto. La sintaxis es una de 

intent(IN) 
intent(OUT) 
intent(INOUT) 

Por ejemplo, considere esta función: 

real function f(x) 
real, intent(IN) :: x 

f = x*x 
end function 
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La intent(IN) especifica que el argumento ficticio (sin puntero) x nunca se puede definir o perder 
su definición a lo largo de la función o su inicialización. Si un argumento ficticio de puntero tiene el 
atributo intent(IN), esto se aplica a su asociación. 

intent(OUT) para un argumento ficticio no puntero significa que el argumento ficticio se con¬ 
vierte en indefinido en la invocación del subprograma (a excepción de cualquier componente de un 
tipo derivado con inicialización predeterminada) y se debe establecer durante la ejecución. El argu¬ 
mento real pasado como argumento ficticio debe ser definible: no se permite pasar una constante 
con nombre o literal, o una expresión. 

De forma similar a antes, si un argumento ficticio de puntero es intent(OUT) el estado de 
asociación del puntero se vuelve indefinido. El argumento real aquí debe ser una variable de puntero. 

intent(INOUT) especifica que el argumento real es definible y es adecuado para pasar y devolver 
datos del procedimiento. 

Einalmente, un argumento ficticio puede estar sin el atributo de intent. Tal argumento ficticio 
tiene su uso limitado por el argumento real pasado. 

Por ejemplo, considere 

integer :: i = 0 
cali sub(i, .TRUE.) 
cali subd, .FALSE.) 

end 

subroutine sub(i, update) 
integer i 

logical, intent(in) :: update 
if (update) i = i+1 
end subroutine 

El argumento i puede tener ningún atributo de intent que permita las dos llamadas de subrutina 
del programa principal. 

9.6. Haciendo referencia a un procedimiento 

Para que una función o subrutina sea útil debe ser referenciada. Se hace referencia a una subru¬ 
tina en una instrucción de cali 

cali sub(. . .) 

y una función dentro de una expresión. A diferencia de muchos otros idiomas, una expresión no forma 
una declaración completa, por lo que una referencia de función se ve a menudo en una instrucción 
de asignación o se usa de alguna otra manera: 


X = fuñe (...) 
y = 1 + 2*func(...) 

Hay tres formas de designar un procedimiento al que se hace referencia: 

1. como el nombre de un procedimiento o puntero de procedimiento 

2. un componente de procedimiento de un objeto de tipo derivado 
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3. un nombre de enlace de procedimiento vinculado 


El primero puede verse como 
procedureO, pointer :: sub_ptr=>sub 

cali subO ! With no argument list the parentheses are optional 

cali sub_ptr() 

end 

subroutine subO 
end subroutine 

y los dos últimos como 

module mod 
type t 

procedure(sub), pointer, nopass :: sub_ptr=>sub 
contains 

procedure, nopass :: sub 
end type 

contains 

subroutine subO 
end subroutine 

end module 

use mod 
type(t) X 

cali x%sub_ptr() ! Procedure component 

cali x%sub() ! Binding ñame 

end 

Para un procedimiento con argumentos ficticios, la referencia requiere argumentos reales correspon¬ 
dientes, aunque es posible que no se proporcionen argumentos ficticios opcionales. 

Considere la subrutina 

subroutine sub(a, b, c) 
integer a, b 
integer, optional :: c 
end subroutine 

Esto puede ser referenciado de las siguientes dos maneras 

cali sub(l, 2, 3) ! Passing to the optional dummy c 

cali sub(l, 2) ! Not passing to the optional dummy c 
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Esto se denomina referencia de posición : los argumentos reales se asocian según la posición en las 
listas de argumentos. Aquí, el dummy a está asociado con 1 , b con 2 y c (cuando se especifica) con 

3 . 

Alternativamente, la referencia de palabras clave se puede usar cuando el procedimiento tiene 
una interfaz explícita disponible 

cali sub(a=l, b=2, c=3) 
cali sub(a=l, b=2) 

que es lo mismo que el anterior. 

Sin embargo, con palabras clave los argumentos reales se pueden ofrecer en cualquier orden 

cali sub(b=2, c=3, a=l) 
cali sub(b=2, a=l) 

Se pueden usar referencias posicionales y de palabras clave, 
cali sub(l, c=3, b=2) 

siempre que se proporcione una palabra clave para cada argumento después de la primera aparición 
de una palabra clave 

cali sub(b=2, 1, 3) ! Not valid: all keywords must be specified 

El valor de la referencia de palabras clave es particularmente pronunciado cuando hay múltiples 
argumentos ficticios opcionales, como se ve a continuación si en la definición de subrutina anterior 
b también eran opcionales 

cali sub(l, c=3) ! Optional b is not passed 

Las listas de argumentos para procedimientos vinculados a tipo o punteros de procedimiento de 
componente con un argumento pasado se consideran por separado. 
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Capítulo 10 


Procedimientos intrínsecos 

10.1. observaciones 

Muchos de los procedimientos intrínsecos disponibles tienen tipos de argumentos en común. Por 
ejemplo: 

■ un argumento lógico MASK que selecciona elementos de las matrices de entrada para ser pro¬ 
cesados 

■ un argumento escalar entero KIND que determina el tipo de resultado de la función 

■ un argumento entero DIM para una función de reducción que controla la dimensión sobre la 
que se realiza la reducción 
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10.2. Uso de PACK para seleccionar elementos qne cumplan una 
condición 

La función de pack intrínseco pack una matriz en un vector, seleccionando elementos basados 
en una máscara dada. La función tiene dos formas. 

PACK(array, mask) 

PACK(array, mask, vector) 

(Es decir, el argumento del vector es opcional). 

En ambos casos, el array es una matriz, y una mask de tipo lógico y compatible con el array 
(ya sea un escalar o una matriz de la misma forma). 

En el primer caso, el resultado es una matriz de rango 1 de tipo y los parámetros de tipo de 
array con el número de elementos siendo el número de elementos verdaderos en la máscara. 

integer, allocatable :: positive_values(:) 
integer :: values(5) = [2, -1, 3, -2, 5] 
positive_values = PACK(valúes, values>0) 

resulta en positive_values siendo la matriz [2,u3,u5]. 

Con el argumento vector rango-1 presente, el resultado ahora es el tamaño del vector (que 
debe tener al menos tantos elementos como valores reales en la mask. 

El efecto con vector es devolver esa matriz con los elementos iniciales de esa matriz sobrescritos 
por los elementos enmascarados de el array. Por ejemplo 


integer, allocatable :: positive_values(:) 

integer :: values(5) = [2, -1, 3, -2, 5] 

positive_values = PACK(values, values>0, [10,20,30,40,50]) 

resulta en [2,3,5,40,50]upositive_values siendo la matriz [2,3,5,40,50]. 

Se debe tener en cuenta que, independientemente de la forma de los argumentos del array, el 
resultado es siempre una matriz de rango 1. 

Además de seleccionar los elementos de una matriz que cumplen una condición de enmascara¬ 
miento, a menudo es útil determinar los índices para los cuales se cumple la condición de enmasca¬ 
ramiento. Este lenguaje común se puede expresar como 

integer, allocatable :: indices(:) 
integer i 

indices = PACK([(i, i=l,5)], [2, -1, 3, -2, 5]>0) 
dando cómo resultado indices son la matriz [1,3,5]. 
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Capítulo 11 

Programación orientada a objetos 


11.1. Definición de tipo derivado 

Fortran 2003 introdujo el soporte para la programación orientada a objetos. Esta característica 
permite aprovechar las modernas técnicas de programación. Los tipos derivados se definen con la 
siguiente forma: 

TYPE [[. attr-list] :: ] ñame [(name-list)] 

[def-stmts] 

[PRIVATE statement or SEQUENCE statement] 

[component-definition] . . . 

[procedure-part] 

END TYPE [ñame] 

dónde 

attr-list una lista de especificaciones de atributos 
ñame el nombre del tipo de datos derivado 

name-list una lista de nombres de parámetros de tipo separados por comas 

def-stmts una o más declaraciones INTEGER de los parámetros de tipo nombrados en la lista de 
nombres 

component-definition una o más declaraciones de declaración de tipo o instrucciones de puntero 
de procedimiento que definen el componente de tipo derivado 

procedure-part una declaración CONTAINS, opcionalmente seguida de una declaración PRIVATE y 
una o más instrucciones vinculantes de procedimiento 

Ejemplo: 

type shape 
integer :: color 
end type shape 
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11.2. Tipo de Procedimientos 


Para obtener un comportamiento similar a una clase, el tipo y los procedimientos relacionados 
(subrutina y funciones) se colocarán en un módulo. 

Ejemplo: 


module MShape 
implicit none 
prívate 

type, public :: Shape 
prívate 

ínteger :: radíus 
contaíns 

procedure :: set => shape_set_radíus 
procedure :: prínt => shape_prínt 
end type Shape 

contaíns 

subroutíne shape_set_radíus(thís, valué) 
class(Shape), íntent(ín out) :: self 
ínteger, íntent(ín) :: valué 

self’/radíus = valué 

end subroutíne shape_set_radíus 

subroutíne shape_prínt(thís) 
class(Shape), íntent(ín) :: self 

prínt *, ’Shape: r = ’, self%radíus 
end subroutíne shape_prínt 
end module MShape 

Más adelante, en un código, podemos usar esta clase de forma de la siguiente manera: 

! declare a variable of type Shape 
type(Shape) :: shape 

! cali the type-bound subroutíne 
cali shape'/oset (10) 
cali shape'/print 


11.3. Tipos derivados abstractos 

Un tipo extensible puede ser abstracto 

type, abstract :: base_type 
end type 
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Tal tipo derivado nunca puede ser instanciado, tal como por 

type(base_type) ti 
allocatable(type(base_type) :: t2) 

pero un objeto polimórfico puede tener esto como su tipo declarado 
class(base_type), llocatable :: ti 
o 

function f(tl) 
class(base_type) ti 
end function 

Los tipos abstractos pueden tener componentes y procedimientos de tipo unido 

type, abstract :: base_type 
integer i 
contains 
procedure fuñe 

procedure(fun_iface), deferred :: def_func 
end type 

El procedimiento def_func es un procedimiento de tipo diferido con interfaz func_iface. Este 
tipo de procedimiento diferido debe ser implementado por cada tipo de extensión. 

11.4. Extensión de tipo 

Un tipo derivado es extensible si no tiene el atributo bind de sequence. Tal tipo puede ser 
extendido por otro tipo. 

module mod 

type base_type 
integer i 

end type base_type 

type, extends(base_type) :: higher_type 
integer j 

end type higher_type 
end module mod 

Una variable polimórfica con el tipo base_type declarado es compatible con el tipo hig-her_type 
y puede tenerlo como tipo dinámico. 

class(base_type), allocatable :: obj 
allocatable(obj, source=higher_type(1,2)) 

La compatibilidad de tipos desciende a través de una cadena de hijos, pero un tipo puede extenderse 
solo a otro tipo. 

Un tipo derivado extendido hereda los procedimientos de tipo enlazado del padre, pero esto 
puede ser anulado. 
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module mod 


type base_type 
contains 

procedure :: sub => sub_base 
end type base_type 

type, extends(base_type) :: higher_type 
contains 

procedure :: sub => sub_higher 
end type higher_type 

contains 

subroutine sub_base(this) 
class(base_type) this 
end subroutine sub_base 

subroutine sub_higher(this) 
class(higher_type) this 
end subroutine sub_higher 

end module mod 

program prog 
use mod 

class(base_type), allocatable :: obj 
obj = base_type() 
cali obj'/oSub 

obj = higher_typeO 
cali obj’/osub 
end program prog 

11.5. Constructor Type 

Se pueden crear constructores personalizados para los tipos derivados utilizando una interfaz 
para sobrecargar el nombre del tipo. De esta manera, los argumentos de palabras clave que no 
corresponden a componentes pueden usarse al construir un objeto de ese tipo. 

module ball_mod 
implicit none 

! only export the derived type, and not any of the 

! constructors themselves 

private 

public :: ball 
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type :: ball_t 
real :: mass 
end type ball_t 

! Writing an interface overloading ’ball_t’ allows us to 

! overload the type constructor 

interface ball_t 

procedure :: new_ball 

end interface ball_t 

contains 

type(ball_t) function new_ball(heavy) 
logical, intent(in) :: heavy 

if (heavy) then 
new_ball%inass = 100 
else 

new_ball°/omass = 1 
end if 

end function new_ball 

end module ball_mod 

program test 
use ball_mod 
implicit none 

type(ball_t) :: football 
type(ball_t) :: boulder 

! sets footballXmass to 4.5 
football = ball_t(4.5) 

! calis ’ball_mod::new_ball’ 
boulder = ball_t(heavy=.true.) 
end program test 

Esto se puede usar para hacer una API más ordenada que usar rutinas de inicialización separadas: 

subroutine make_heavy_ball(ball) 
type(ball_t), intent(inout) :: ball 
ball'/omass = 100 

end subroutine make_heavy_ball 


cali make_heavy_ball(boulder) 
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